A new quantum mechanical photon counting distribution formula 
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| Photon counting by photoelectron devices is a very important topic in quantum optics. Since it can judge the 
. nonclassical features of light fields, and most measurements of the electromagnetic field are based on the absorption 
' of photons via the photoelectric effect [H, 0). So the problem of photoelectric detection attracts the attention of 
t-H , many physicists and scientists. Expression of quantum mechanical photon counting distribution was first derived by 
Kelley and Kleiner [H . For the case of single radiation mode the probability distribution p (m,T) of registering m 
photoelectrons in the time interval T is given by0-Q 



p (m, T) = Tr 



(1) 



where p is a single-mode density operator of the light field concerned, : : stands for normal ordering, and £ is called 
the quantum efficiency (a measure) of the detector. Mogilevtsev's group has as well discussed no counts registered 
on a photonic detector Q, i.e., for the case of to = in Eq.([T]). Recently, Fan and Hu Q derived two formulas for 
photocounting; one involves p's coherent state mean value and the other is related to p's Wigner function. 

As an example, when p— \a) (a\ is a pure coherent state, \a) — exp (cva^ — a* a) |0) = D (a) |0) , and a\a) — a \a), 



then the probability distribution is 



HrKyU-*^: \a)=W°LLl . (2) 



which is the Poisson distribution as expected. Due to the operator identity e A ° ta =: exp[(e A — f) a* a] : Q and 

e ataIn(l-e) o m e -otaln(l-0 = a m (1 _ £)" m , (3) 



as well as 



jm a m = N (N —!)•■• (N — m+ 1) , N = a f a. (4) 



we see 
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= £ Lk m (i-r m i«>H. (5) 
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where |n) = ^— |0) is the number state. Substituting Eq.© into Eq.([T]) yields 

p(m,T)= E^»U W-0"~ m , (6) 

n—ra ^ ' 

where "P„ = (n| p |n) . Eq.Q is named as Bernoulli distribution in particle-number representation [4, 5]. However, for 
many density operators other than the pure coherent state, directly using Eq.([T|) is difficult and is inconvenient, so in 
this work we shall reform it. 

II. NEW PHOTON COUNTING DISTRIBUTION FORMULA DERIVED BY DENSITY OPERATOR'S 

P- REPRESENTATION 

Now we instead use density operator's P-representation0 

p= J d 2 aP{a) \a) (a| (7) 

to analyse Eq. ([!]), \a) is the coherent state. Substituting Eq.([7]) into Eq.([T]) yields 



p(m,T) = / d z aP(a){a\: 



£a T a , 
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= /dV^MT^H 2 . (8) 



Using P (a)'s expression [h 



P ( a ) = e H 2 / *l (~(3\ p |/3) e \P\ 2 +^'-^' , ( 9 ) 

7T 



where |/3) is also a coherent state, Eq. ([SJ becomes 



m! / 7r 



(£-1) 

where L m (x) is the Laguerre polynomial, 



<-«<■!»• (J^). w 



i=0 1 ' i'=0 

and we have used 



— ^z"z* m e e ' z ' = e~ BC l A Y ° - .- T ,Re6<0. (12) 

77 n(n-0!(m-Z)!(-e) n+m -' +1 



Eq. (fT0|) is a new formula, different from that in Ref.Q- 



III. APPLICATIONS 

As some applications of Eq. ([T0l) . substituting the pure coherent state p= \a) (a\ into Eq. fTO)) . we have 



3 



where (a \ j3) = exp [— | (|a| 2 + \f3\ 2 ) + a*0] ■ Enlightened by the generating function of the Laguerre polynomial [TTJj 



so 



1 d n 

L n (x) = — — — 
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we rewrite Eg. (1131) as 



P|«><a| (m,T) 
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where we have used the formula 
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(17) 



It is clear that Eq.(fT6|) exactly agrees with Eq.(|2]). Comparing Eqs. (fl"3|) with Eq. ([T6|) implies a new integration formula 



7T \4 — 1 

Especially, when a — 0, in Eq. ([18[) only m = term survives, so we obtain 

|^e^X(|/3| 2 ) = l. 

Now we calculate another density operator 



p c = (1 - e-f) 



with / = -^y, which represents a filtered one- mode chaotic light. Due to 

(1 - e-f) (-P\ : e( e_/ - 1 ) ata : |/3) = (l - e~f) e< l+e ~ , )^\ 
substituting Eq. (|2"Tj) into Eq. fTU)) . we have 

«|2 



(18) 
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(21) 
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With the help of the formula [IJ] 



drL m (r) e- 9r = 



(g-i) ? 



Equation f|22|) becomes 



Pc (m,T) 
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m+l ' 



(22) 



(23) 



(24) 
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Noting (e f - 1) 1 = (e^r — 1^ = n, representing the Bose-Einstein statistics, Eq. ([24"l) leads to the photon counting 
probability for the chaotic light 



pc(m,T) 



m+l ' 



which agrees with the known result, so our approach's correctness is confirmed. 
When p s is a pure squeezed state 

tanh A „ 1 2 , . , , tanh A „ 2 

p s = sechAe^- a |0) (0| e~^ a , 
where A is the squeezing parameter, after p s being substituted into Eq.(9) we see 



Ps (m,T) 



\m+l 



Then using Eq. ([T5j) and the integral formula[l 
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with Rc(C + / + g) < 0, Re(g^) < or Re(C - / - g) < 0, Re(gf^f ) < 0, we have 



(25) 



(26) 



(27) 



(28) 



Ps (m,T) 
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where 



G 2 = (£ - l) 2 tanh 2 A. 
Using the generating function of Legendre polynomial (Tl| 



(29) 
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finally we have 
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which relates to the Legendre function, and this is a new result. 
For the last example, when p^ is the displaced chaotic field 



p d = (1- e~ f ) D (a) e-t'^D- 1 (a) 



(33) 



(34) 
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with D (a) = exp (aa) — a*a), we convert pd into its normal ordering formflij 



Pd 



n + 1 



exp 
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n+1 



(35) 



and substitute it into Eq. (fTU|) yields 



<Pd(m,T) 
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(36) 



where we have considered Eqs. (fT5|) and (ITTl) . It is interesting to see that when £ = l,n = —1/2, /?d in Eq. ([35|) becomes 
Pd = 2: exp [—2 (a — a) (a* - a*)] : = 27rA(a), where A(a) is the Wigner operator. Then using : e~ a a : = |0) (0| , 
|m) = ^7= 10), we have 

1 ' \/m 1 ' 



p d (m, T) = 2tt (m| A(a) |m> = 2 (-1)™ e- 2 H 2 L r , 



4|a| 



2ttW (a) , 



where (a) is just the Wigner function of the number state |m)[l5|. This is a good check for the correctness of 
Eq.([56]l. 

In summary, by virtue of density operator's P-representation in the coherent state representation, we derive a new 
quantum mechanical photon counting distribution formula. As its application, we find that the photon counting 
distribution for the pure squeezed state relates to the Legendre function. 
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